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Abstract 

Within the NRQCD factorization approach the production of heavy 
quarkonia can be calculated in perturbative QCD in terms of a few 
production-process-independent NRQCD matrix elements. We test 
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ing their values as determined from charmonium production at pp, 
ep, and e + e~ colliders, and in Z°, T, and B decays. 
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1 Introduction 



The production of heavy quarkonia in high-energy reactions can be calculated 
in perturbative QCD in terms of a few non-perturbative parameters that are 
independent of the production process. More precisely, a factorization theorem 
holds, which expresses the cross section to produce a quarkonium H as a series 
of terms 

a(H)=J2cr(QQ[n]) (0 H [n]) , (1) 

n 

where a(QQ[n]) is the cross section to produce a quark-ant iquark pair with small 
relative momentum at distances 1/m or smaller (m is the heavy-quark mass). 
The quantity (O n [n]) gives the probability for the heavy-quark pair of state n to 
bind into the meson H. On the one hand, the binding occurs at distances much 
larger than 1/m, ensuring the existence of the factorized form (p. On the other 
hand, the binding energy for neither charmonia nor bottomonia is large compared 
with the dynamical low-energy QCD scale Aqcd- Hence the long-distance matrix 
elements (MEs) (0 H [n\) cannot be computed perturbatively. However, there 
exists a power-counting scheme that provides a hierarchy of the MEs so that to 
any desired accuracy only a given number of terms contribute in (|I]). 

The factorization ([]]) has been derived perturbatively using the so-called non- 
relativistic QCD (NRQCD), an effective field theory appropriate for mesons con- 
taining a heavy-quark pair It therefore holds to leading order in Aqc^/m. 
Higher-twist corrections occur both in the production of the heavy-quark pairs 
and in their subsequent fragmentation into the meson. The importance of the 
NRQCD MEs (0 H [n\) is governed by their scaling with v — \v\, the relative ve- 
locity between the heavy quark Q and the heavy antiquark Q within the bound 
state. Equation [I] is hence a double expansion in v 2 and a s (= a s (/x) with jj, ~ m). 
The occurrence of v 2 as an expansion parameter is in line with the success of 
non-relativistic potential models in describing quarkonium spectroscopy and dis- 
tinguishes quarkonia from heavy-light systems where Aqcd/^ is the expansion 
parameter. 

In the non-relativistic limit the quarkonium is a pure QQ state with quantum 
numbers that match those of the meson, for example a cci( 3 <Si) state in the case of 
the J/ip (J PC = 1 ) and a cc"i( 3 Pj) state in the case of the J ++ Xcj mesons. The 
phenomenologically most important consequences of the NRQCD factorization 
approach (FA) are two-fold. First, two MEs contribute to Xcj production at 
leading order in v 2 , namely (0* cJ ( 3 P/)) and (0$ cj ( 3 Si)) . That is to say, besides 
the contribution from the leading Fock state, the colour-singlet cc( 3 Pj) state, 
there is, at leading order in v 2 , a contribution where the heavy-quark pair is 
produced in a colour-octet state. 

The second, perhaps more dramatic, consequence concerns the 1 mesons 
(J/ip, ip', T(raS)): although there is just one contribution corresponding to the 
short- distance production of the leading QQi( 3 S\) Fock state for v 2 — > 0, rel- 
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ativistic corrections are sizeable, even partly dominating the production cross 
section, at least in the charm system: v 2 ~ 0.3 is not small enough for gluon 
radiation from an almost point-like QQ pair produced in the hard interaction 
(governed by a 8 ) to dominate gluon radiation from the pair at later times, when 
it has already expanded to the quarkonium size ~ l/(mv ) (governed by v 2 ); com- 
pare, for example, (^) and (||) below. Therefore contributions to the production 
of 1 mesons are important where the heavy-quark pair is produced at short 
distances in a colour-octet state. Colour-octet MEs first appear at relative order 
f 4 . Different quarkonium-production reactions differ in the relative weights in 
which these MEs enter the 1 cross sections, i.e. by their short-distance parts 
a(QQ[n}) with n = 3 S[ 1] , 3 S{ 8 \ 1 S^ S \ and 3 Pj 8) up to and including 0(v 4 ). 



2 Charmonium production at the Tevatron 

In pp collisions at high energies, charmonium particles come from the decay 
of b-flavoured hadrons and prompt production, the latter consisting of direct 
production and feed-down from higher charmonium states (except for if)', the 
highest-lying charmonium state). Short-distance processes that produce direct 
ip(nSys in the leading colour-singlet state are 

g + g -> ccxfSO+g : a s 3 ^ , (2) 

g + g -> g + g", g*->cci( 3 Si)+g + g : al\ , (3) 

P± 

where the power of a s and the p± fall-off (at 90° in the partonic c.m.s.) have 
been indicated. At large p±, (|J) dominates and can be calculated by folding the 
gg — > gX subprocess with the g — > ipX fragmentation function 0. With values 
of the colour- singlet MEs as obtained from potential model calculations, (0) and 
(|3[) underestimate the observed yield of direct J/ip and mesons by a factor 
~ 30. The data can, however, be accounted for (Fig. [I]) if the short-distance 
production of colour-octet QQ states is included [§, [§] 

Attl 2 

g + g -> cc 8 [n]+g n = 1 S , 3 P J : a 3 s v 4 — —, (4) 

PI 

g + g - g + g", g^ccg^) : alv'\. (5) 

P± 

The p± shapes of these two mechanisms are sufficiently different over the mea- 
sured range to isolate the 3 Si colour-octet ME, but the remaining two MEs can 
be extracted only in the linear combination 

Mf = (OfC^o)) + ^ K( 3 ^o)) , (6) 
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Figure 1: Fits to direct J/ip production at the Tevatron ||. 

where a pa 3.5. The numerical values are indeed consistent with the expectation 
from velocity scaling, see table [I]. The values are, however, subjected to a number 
of uncertainties: 

• The value of the charm-quark mass; it affects mainly the overall normaliza- 
tion and should always be supplied when quoting values for the ME. 

• The parametrization of the parton-distribution functions (PDF) and the 
choice of the factorization scale; both influence the normalization as well as 
the shape of the p± distribution. 

• The value of Aqcd (which is not independent of the PDF) and the renor- 
malization scale fiR, note that processes (f|) and (||) start at high order in 
a s . The values of Aq CD and hr affect also the p± shape, since fi 2 R ~ m 2 +p]_. 

• The lack of higher-order perturbative QCD corrections. Particularly im- 
portant are: (i) initial-state radiation and intrinsic kj_ effects since they 
may modify the slopes of (f|) and (^) considerably M; (ii) Colour-singlet 
production via the 0(at) reaction g + g — > cc"i( 3 Si) +g + g, since it has the 
same p ± fall-off as (gj) [0 . 

• Higher-order v 2 corrections are important close to boundaries of partonic 
thresholds ||. 
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Table 1: Scaling of NRQCD matrix elements with v (modulo v 3 ) and numerical 
values in GeV -3 as obtained from a potential-model calculation (second and third 
row) || and fits to the Tevatron p± distribution (last two rows) || . 

• Higher-twist corrections. At large p± these are probably small for the short- 
distance cross sections. However, sizeable corrections may be present for 
the fragmentation functions at the input scale. 

• Values of the colour-singlet matrix elements; at the Tevatron this is impor- 
tant for Xcj production only, i.e. the extracted value of (0^ CJ ( 3 Si)). 

From the above discussion we expect the current determination (table |T|) of 
(Og^ nS \ 3 S\)) to be on the low side, Mf^ 1 ^ to be overestimated, and a strong 
correlation between (0f cl ( 3 Si)) and {Of 1 ( 3 P 1 )). 

Although it is reassuring that the FA can accommodate the data, the result 
is not too surprising given that we have two additional mechanisms that scale 
as l/p± and l/pj_, respectively, with free normalizations. Additional consistency 
checks are therefore needed. The polarization of the J/ip, measurable via its 
decay into is one such test since it is uniquely predicted without the need of 

additional long-distance MEs up to 0(v 4 ) |4]]. Polarization measurements should 
soon become available at the Tevatron. Here we focus on another aspect of the 
FA, namely the universality of the ME. To this end we compare MEs extracted 
from various reactions. 

3 Z-decay 

The Feynman diagrams for quarkonium production in e + e~ annihilation are 
shown in Fig. |2|. Results of a recent prediction |I| are compared with LEP 
data ]11], |12| in table |^. At e + e~ cm. energies i/i that are large with respect 
to the heavy-quark mass, for example in Z° decays, contributions with addi- 
tional gluons (lower diagrams in Fig. ||) are suppressed by powers of mj y/s. The 
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Figure 2: Feynman diagrams for quarkonium production in Z° decays. 



dominant contributions arise from QQ production in the leading (colour- singlet) 
Fock state (upper left diagram in Fig. |2|) and the production of a QQs( 3 Si) pair 
through gluon splitting (upper right diagram in Fig. 0). The former contribution 
is the only one present in the so-called colour-singlet model (CSM). As table ^| 
reveals, this colour- singlet contribution alone seems low with respect to the J/V> 
data, although still compatible, within errors. 

The data are well described, after inclusion of the colour-octet part, with a 
value for (Og^( 3 Si)) as extracted from fits to the Tevatron p± distribution. In 
the case of T production, the statistics per experiment is too low to allow for 
definite conclusions. Further improvement can be expected because not all Z° 
data have yet been analysed. Since the error of the data is statistics-dominated, 
a combined LEP analysis would be welcome. Information on colour-octet MEs 
other than (0^( 3 Si)) will come only from e + e annihilation at lower energies, 
BES, CESR, and B-factories. 

The LEP data are also consistent with the expectation of the colour- 
evaporation model (CEM), a model of quarkonium formation that assigns pre- 
cisely one process- independent, long-distance factor f H = (0 H [n\) to each 
quarkonium particle H |TJ[. The cross section a(QQ[n\), see ([I]), is obtained 
by integrating the QQ production cross section up to the threshold of the pro- 
duction of a pair of heavy-light mesons. With parameters f H extracted from 
fixed-target hadroproduction, the CEM is also able to describe the Tevatron p± 
distributions . 
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CSM CEM NRQCD 


DELPHI OPAL ALEPH L3 


J/lf; 


7.8 23 26 


44±|g 19 ± 10 30 ±9 27 ±12 




1.7 1.7 6.6 


< 124 10 ±5 < 7.3 < 7.6 



Table 2: Branching ratios and 95% CL upper limits in 10 -5 of Z° decays: pro- 
duction of the sum of the three lowest-lying T states and prompt J /if) production 
(i.e. excluding B decays but including feed-down from Xcj and if)'). 

4 Low-energy electron— positron annihilation 

In low-energy e + e~ annihilation, J /if) production is unique in the sense that it is 
the only process known so far where the colour-octet 1 5'o and 3 Pq MEs enter with 
widely different weights. At low energies, the dominant contributions involve a 
cc pair plus one or two perturbative gluons (lower Feynman diagrams in Fig. |2|) 
l|, |l6fl . The cross sections are readily evaluated and results presented in table ^| 



and Fig. |3|. 

Observe first that the ME combination @ enters with a varying between 11.3 
and 3.8 in 5 < ^/s < 10.6 GeV. The central values of the experimental data yield 
for a s = 0.278, m c = 1.5 GeV, and (£>^( 3 Si)) = 1.16 GeV 3 

(o^'CSo)) = 4.0 x 10- 4 GeV 3 , V 8 , ° V = 9.3 x lO" 3 GeV 3 , (7) 

and hence = 3.3 x 10~ 2 GeV 3 (unprimed curves in Fig. |3|). The last number 

is somewhat smaller than the naive Tevatron estimate, see table [I]. 



Our numbers are similar to those found in [[16] . There it was claimed that 
(Og^( 3 P ))/m 2 can be extracted very precisely from low-energy e + e~ data while 
(Og ( 1 S'q)) is hardly constrained, based on the dominance of i? 8 ( 3 -Po) over R\ 
and R^Sq) at low energies and the fact that i? 8 ( 1 S , ) is a small contribution for 
all ^s. Indeed, varying a s in (0.24, 0.30), m c in (1.4, 1.6) GeV, and (C 1 J/v '( 3 5i)) in 
(0.80, 1.52) GeV 3 yields, for the central experimental values, the following rangesQ: 

- 1-2 < 7^3 (O^CSo)) < 1.2 , 1.1 > — ^ (oi /lp ( 3 P )) > 0.78 , (8) 



GeV \ v 'I " GeV 3 m? 
implying 2.6 < 10 2 M J 3 $ / 'GeV 3 < 4.5. 



1 These are the maximal ranges obtained by varying one parameter at a time; the uncertainty 
associated with each of the three parameters is about the same. 
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0.2 


56 


20 


3.8 


39 


-23 


10.6 GeV 


10 ±4 


3.4 (Ol /lp ( 3 S 1 )) 


170 (Ol'+fSo)) 




3.9 


0.07 


6.0 


14 


3.9 


12.5 


-2.5 



Table 3: Cross section of prompt J/ip production in e + e~ annihilation (excluding 
b decays) in units of lOV^: Data [O (i?[exp.]), colour- singlet contribution (R\), 
and colour-octet contributions through 1 S and 3 Pj intermediate cc pairs. The 
PLUTO data a = (31 ± 21) pb have been multiplied by (1 - 0.23) 6/7 to remove 
the if)' contribution and to update the leptonic branching ratio J/ip — > [i + . 
The theoretical values are for a s = 0.278, m c = 1.5 GeV, the colour-singlet ME 
as given in table [I] and the colour-octet MEs (0). Also given are the results of 
one of the extreme solutions © . 

However, a few remarks are in order. First, the large experimental er- 
rors forbid definite conclusions^. Even for fixed theoretical parameters a s , m c , 
(Pi( 3 Si)) only the following, much wider range is obtained: 

- °- 73 < ^3 (O^CSo)) < 7.4 , 2.2 > (O^CP )) > -0.38 , 

(9) 

implying 

M 3 J.t = (3.3 ± 2.7) x 10" 2 GeV 3 . (10) 

Insisting in (0%^ S °^) > 0, the minimally allowed value of M. 3 /^ increases from 
0.57 x lO- 2 GeV 3 to 1.1 x 10" 2 GeV 3 . One pair of MEs fulfilling (g) ((£>g Ml5o) ) = 
0.05 GeV 3 , (Ol^( z P ))/m 2 c = 9.5/10 5 GeV 3 ) is shown in Fig. | (primed curves). 
If we allow in addition the (independent) variation of the theoretical parameters, 
even negative values of are not excluded. 

Negative values for (O^ ( 3 Pq)) are not per se excluded, although at first 
sight these would lead to negative cross sections at low cm. energies since the 
cross section of ccg( 3 Pj) production grows at low energies while all others vanish. 

2 Recall that PLUTO f7) observes an excess of only 4 prompt c + c — ► J/ipX events; out of 
these, only the two events with both decay muons reconstructed have been used here. 
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Figure 3: Fits to direct J/?/> production in e + e annihilation. 



However, the e + e~ — > J/ipX cross section at low energies is no longer given by 
the fixed-order calculation. The reason is the singular behaviour of the 3 Pj cross 
sections: for £ = 1m c j ^fs — > 1, these diverge as 1/(1 — £). Clearly this indicates 
the breakdown of the fixed-order in v 2 calculation: t> 4 /(l — £) is no longer a small 
parameter and higher-order terms in v 2 grow more rapidly as £ — > 1. The energy 
dependence at low energies can only be described once these terms of higher order 
in v 2 are resummed. Eventually, as 1 — £ < v 2 , higher-twist contributions become 
important, too. Of course, the result of the fixed-order calculation can still be 
used if we smear over a sufficiently wide region in ^/s. 
Near threshold the cross section is given by 



e + e 



« Uo^Cs,) 



direct 



a s 

+ — 

7T 



,(i-O + G)'f?<i-0 



i i 



"c J=0 



, (11) 



where do = 7i 3 a 2 m el/(6m^), e c = 2/3. The apparent singularity at £ = 1 can be 
absorbed through the factorization-scale dependence of the leading colour-singlet 
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ME owing to the evolution equation derived in JI], [T8 



£ o 8 ( 3 /',)-^Pi( 3 s 1 i . 



c j=o 



37rm^ 



(12) 



Hence a finite cross section is arrived at when averaging the threshold region over 
a range of order m c . 



5 B decay 

The B meson is known to decay into charmonia with branching ratios between 0.1 
and 1% |19 |. Branching ratios into charmonia of the other b-flavoured hadrons 
are not known. The hard process is given by 

b — > cc c [n] + X 8 , n = S , 3 Si , 3 Pi , D— waves , . . . , (13) 

where X s is a system containing a strange quark and, as usual, the subscript 
c indicates a colour-singlet (c = 1) or colour-octet (c = 8) heavy-quark pair. 
Calculations of B decays into charmonium states based on ( |T3[ ) are valid to leading 
order in v 2 and leading order in Aqco/^b- Using the results of [pH[] , we can 
find the B into charmonium decay-branching ratios in terms of the semileptonic 
branching ratio 

9 mi f(m c /m h ) 
x{[2C + -C_] 2 ff + [CV + C_] 2 ^ff} . (14) 

Here f(r) is the phase-space function for the b semileptonic decay 

f{r) = l-8r 2 + 8r 6 -r 8 -24r 4 ln(r) , (15) 
R(z) is the phase-space function for b — > H + X, 

R(z) = - \l-(2z) 2 } 2 , (16) 

and 

rf = (i + 8^)(of( 3 ^))--^(pf(X 

x ' 6 mi N 



+ 3(cV'( J 5„)) (l + Sc 3 )^(Of( 8 P0) • (17) 



c 



Of course, for a particular charmonium particle H and a specified colour state c 
not all MEs in (|17|) contribute to any given order in w 2 , see table |l|. The second 
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ME in (f7|), (V?{ 3 Si)), is a v 2 correction to (Of ( 3 Si)) and will be neglected in 
the following. 

Finally, C± are factorization-scale-dependent Wilson coefficients describing 
the evolution from the W-mass scale down to a scale of the order of the b- 
quark mass. Note that the coefficient of the colour-octet contribution is strongly 
enhanced compared with the colour- singlet one: 



3 
2 



' c+ + c_ ' 

2C+-C- 



50 . 



(18) 



In the following we take as central values C + + C- = 2.2, 2C + — C- = 0.38, 
m c = 1.48 GeV, m b = 5.0 GeV, Br [b -> c£ut] = 10.3%. In contrast to the 
colour-octet Wilson coefficient, the colour-singlet one is strongly scale sensitive, 
since it is given as the difference of two large numbers. In order to bracket its 
uncertainty we shall vary 2C + — C- between 0.19 = 0.38/2 and 0.645. With the 
latter value, Br [B — > J/ipX] is saturated with only the colour-singlet contribution 
(see below). Varying the colour-singlet coefficient in (0.19,0.645) gives the range 
from 17 to 201 for the ratio ([18]). 

As is true for most quarkonium production reactions, B decays into charmonia 
are sensitive to only the combination @ of colour-octet ME. Moreover, the value 
of a is similar to the one at the Tevatron, a = 2 (1+8 z 2 ) rs 3.4. This follows imme- 
diately upon using the approximate symmetry relation (Of ( 3 Pi)) ~ 3 (Og ( 3 Po))- 

Restricting to the lowest order in v 2 terms, separately for the colour-singlet 
and colour-octet contributions, we find the following direct branching ratios 



GeV 3 



0.24 -r^w 



Br [B -> V + X] = 
J/V> : 0.80 ± 0.08 0.28 (0.07,0.80) 
if)': 0.34 ±0.05 0.18 (0.04,0.52) 

0.10 



+ 11 GeV a 
0.13 

0.046 
0.024 



+ 21^4 

0.92 

0.38 

0.20 (k = 



(19) 



0.53) 



and 



Br [B^XcJ + X]-- 
Xd ■ 0.37 ±0.07 

Xc2 ■ 0.23 ±0.10 



0.48 



0.070 (0.017,0.20) 
0.046 





+ 12 



GeV 
0.28 

0.19 
0.47 
0.29 



(k = 0.66) 
{k = 0.63) 



(20) 

The experimental numbers (first rows) are the CLEO measurements |H| of the 
direct branching ratios (obtained by subtracting the feed-down from ip' — > x c j X 



10 



and XcJ — ► J/V'T)- The numbers in brackets denote the range of the colour- 
singlet contribution obtained by varying 2C + — C_ in the range specified above. 
The numbers quoted for k 7^ 1 are obtained by using meson masses rather 
than quark masses in the phase-space function, i.e. by multiplying ([[4]) by 
k = R([rriH /2]/ m B) / ' R{m c / m\>) . Clearly there are other v 2 , AQco/^b, and even 
Aqcd/^c corrections besides this trivial phase-space factor. Hence these num- 
bers should be regarded as an only very naive indication of the uncertainties due 
to these corrections. 

We observe that, by stretching the poorly determined colour-singlet Wilson 
coefficient 2C + — C_, production of J /if) and if)' can be explained merely by the 
colour-singlet contribution, in contrast to Xcj production. This is in accord with 
the fact that there are two Xcj MEs in leading v 2 while only one if> ME survives 
for v — > 0. On the other hand, the values of the if> 3 Si colour-octet MEs are well 
compatible with the CLEO data for almost the whole range of 2C + — C_. In 
any case, the colour-octet ME combination A4^ A is too large, for example, by a 
factor of about 2 for the central value of 2C + — C_. This confirms the expectation 
raised above that it has been overestimated in fits to the Tevatron data. 

An inspection of the Xc2 width reveals that the Xcj 3 S\ colour-octet ME is a 
factor of about 2 too large. The Xci width then requires the colour-singlet ME 
(C^ cl ( 3 P 1 )) to be larger, by a factor of about 3 for the central value of 2C + — C_. 
This may well be consistent with the Tevatron data where there is a similar strong 
correlation between the two MEs. 

While the leading-order Xcj ME is rather well known, MEs of higher orders in 
v 2 have not yet been determined. However, relativistic corrections are likely to be 
important for Xcj production in fixed-target experiments and in e + e~ collisions 
^j] . Corrections to 0(v 2 ) arise merely from operators with two extra spatial 



derivatives, but the same quantum numbers as the leading ones. Contributions 
that scale as v 4 relative to (|20|) give the correction 



ABr[B^ XcJ + X] = {0.24 (OffSt)) +0.42 (Of (21) 



where Co is an as yet uncalculated coefficient. Nonetheless, it is clear that B 
decays into Xcj states provide an upper limit on M^i- Applying the scaling 
f 4 ~ 0.1 to the two leading-order MEs, we find values for -M^ 1 varying between 
2 x 10~ 3 and 1.5 x 10 -2 . B-decay data certainly favour values at the lower end. 

Let us finally mention that the Xc2-to-x c i ratio is a serious problem for the 
CEM, which predicts the ratio 5/3. 
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6 Upsilon decay 



Inclusive charmonium production in T decays is another place to test the NRQCD 
factorization approach. In principle, colour-octet contributions in both the bb de- 
cay and the cc production should be taken into account. However, relativistic 
corrections from the bb in a colour-octet state turn out to be small [22|, in ac- 
cordance with the naive expectation of a suppression by (fb/f c ) 4 ~ 1/10 relative 
to a contribution having the cc in a colour-octet state. 

A full calculation of T decays into charmonia does not yet exist. The colour- 
singlet contribution to direct J/ip production starts at 0(a®) and consists of 
the tree diagrams bbi( 3 Si) — > cci( 3 Si) + 4g and the loop diagrams bbi( 3 5i) — > 
cc"i( 3 S'i) + 2g containing box diagrams. Only an estimate of the 0(a^) colour- 
singlet contribution to indirect J/ip production through Xcj decays exists, yielding 
about 1/20 of the measured branching ratio f2~2| , |23f . 

Colour-octet contributions to direct J/^ production in T decays start at 0(a 4 ) 
and are hence enhanced by 1/a 2 . Although calculations of the colour-singlet 
contributions are still lacking, a comparison with data can be used to set upper 
limits on certain colour-octet MEs. The largest colour-octet contribution found so 
far is bb 1 ( 3 S'i) — * ggg* followed by g* — > cc 8 ( 3 Si) p2| . Another potentially large 
contribution is the loop process bbx( 3 Sx) — > cc 8 ( 1 S , , 3 Pj) + g proceeding through 
virtual gluons. This process can be related^ to the radiative decay QQi( 3 Si) — > 
7 + ?9i( 25+1 -^j) calculated some time ago |2~1. We find 



Br [T -> J/-0 + X] = BRx + R 



7r 2 — 9) m 3 
|o.571 (O^CS^ + QM^mS} 
BRx + 1.76 x 10- 2 M / *( 3 S 1 )\ + 0.247 x 10- 2 A^ 



8 V U D / T i ^ iu - ,vt 3.8 

4 



(ll±4)xl0~ 4 = BRx + {1.94 + 1.09} x 10 , (22) 

where BRx denotes the missing part from unknown colour- singlet contributions 
and indirect J/ip production. The second equation has been obtained with R = 
Br [T — > 3g] ~ 0.83, a s = 0.278, m c = 1.5. The number on the LHS of the last line 
is the CLEO measurement [^] and the values of the RHS follow from insertion 
of the ME values in table [l]. The latter are consistent with the measurement if, 
as expected, the colour- singlet part is small and the indirect contributions does 
not exceed, say, 50% of the total J/ip yield. More definite conclusions can only 
be drawn once the missing pieces have been calculated and/or more experimental 
information is available, for example, on the J/ip energy distribution or on the 
Xcj production rate. 



3 Whilc finishing this paper we learnt that this process had just been taken up also in 
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Figure 4: Photoproduction of J/ip compared with HERA[^] data. Theoretical 
curves at ^ = 100 GeV for {o{'^ Sl) ) = 1.16 GeV 3 , (0 J 8 /f{1So) ) = 0.03 GeV 3 , 

(C^( 3 P )) = 9/10 3 GeV 3 , a s = 0.278, and m c = 1.5 GeV, using the leading- 
order GRVJ30] parton distributions with \i = 2m c . The colour-singlet contribu- 
tion is shown also for a s = 0.3, m c = 1.4 GeV (dotted line). 



7 Photoproduction 

Photoproduction of J/ip has been claimed to be in conflict with NRQCD factor- 
ization ||27|| . In fact, the short-distance production of merely the leading colour- 
singlet state via 

T + g-ce^SO + g (23) 

yields a shape of the differential cross section da/dz, in good agreement with the 
§ (Fig. 



HERA data |28[ (Fig. |j), and also in agreement with low-energy data [2S]. Here 
z is the fraction of the photon energy carried by the J/ip in the proton rest frame. 

The absolute normalization suffers from uncertainties similar to the ones listed 
in section |2| for J/ijj production at the Tevatron: valueQ of m c , parametrization 
of the parton-distribution functions, value of Aqcd, choice of factorization and 
renormalization scales^, and initial-state transverse-momentum effects, i.e. effects 

4 Note that the m c dependence is partly compensated by corresponding changes in the ME; 
this has not been taken into account in Fig. |]. 

5 The scale dependence has been reduced through the inclusion of next-to-leading order 
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arising from the intrinsic transverse momentum k± of the gluon within the proton 
and initial-state gluon radiation. The latter effects mainly influence the J/ip 
p_i_-distribution at low p±, but do affect the z distribution as well, mostly its 
normalizations, less so its shape. This occurs in particular if a cut on p± is 
applied as low as 1 GeV, not much larger than (k±) ~ 500 MeV. Recall that 
fixed-order calculations are reliable only at large p± or for p_i_-integrated cross 
sections but do not adequately describe the shape at low p±. 
Relativistic corrections of order v 4 arise from 

7 + g -+ cc 8 [n] + g n = 3 S U 3 Pj 

7 + q -> cc 8 [n]+q n = % 3 St, 3 Pj , (24) 

and from reactions involving the partonic content of the photon. Although such 
"resolved-photon processes" become more important as the cm. energy increases, 
they still affect the small-z region, only, say, below 0.3 at HERA energies. The 
quark-initiated reactions in (j24"D are small with respect to gluon-initiated ones at 
HERA energies. Moreover, 7g — > ccs( 3 5'i) + g is identical to (^3|), but down in 
magnitude by a factor of about 50. Most important are therefore the other two 
gluon-initiated reactions in (|24]). In fact, their contributions seem to be in clear 
conflict with the HERA data when using the MEs as given in table [l], see Fig. f|. 
This is not an artefact of the particular values of the S- and P-wave MEs chosen 
in Fig. d, since both contributions are rather similar in shape. (To be precise, 
the value of a in the ME combination (^) is a = 4.7 ± 1.4 in the range < z < 1 
for p±> 1 GeV.) 

Three reasons make us believe that we do, in fact, not encounter a breakdown 
of NRQCD factorization, namely (i) the size of the ME, the treatment of (ii) the 
hard process and (iii) the hadronization. First, we have seen that all but the 
Tevatron estimates prefer smaller values of the ME combination A4 J J^ 5 . And we 
have given arguments why the value needed to explain the Tevatron data is likely 
to be smaller. 

Second, (f24j) contains contributions that correspond to the evolution of the 
gluon (or quark) distribution functions. These terms arise from diagrams with 
t-channel gluon exchange and have to be removed via mass-factorization in calcu- 



lations of fully p_i_-integrated cross sections [ft^l . Cross sections integrated above 
some p_i_-cut are reliably described by the fixed-order (in a s ) results based on 
(|24|) only if the cut is sufficiently large. If the p± cut is as low as, say, 1 GeV 
then the average p± of the fixed-order calculation is considerably smaller than 
what would be obtained from a proper calculation including soft-gluon (initial- 
state) radiation. Since contributions from small p± lead, on average, to harder z 
distributions and vice versa, we expect a flatter z distribution after inclusion of 
initial-state radiation and intrinsic k i . 



perturbative corrections [31|. 
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The third aspect concerns the sensitivity of the z distribution to energy ranges 
in the ccg[n] — > J/ip transition of the order of the binding energy (and smaller). 
The physical reason is clear: the momentum carried away by light hadrons in 
the hadronization process makes it very unlikely that z values close to 1 can be 
reached. Qualitatively, the hadronization causes a smearing of the J/tp momen- 
tum over a region 5z ~ v 2 ~ 0.25-0.30, 5p± ~ mv 2 ~ 0.5 GeV. Technically, 
the NRQCD expansion parameter at z close to one isQ v 2 /(l — z) rather than 
v 2 ||. The amount of smearing caused by this breakdown of the fixed-order (in 
v 2 ) calculation makes the latter unsuitable to predict the z distribution. Predic- 
tivity can, however, be restored upon introduction of universal shape functions. 
Attempts computing these shape functions in the Coulomb limit (m — > oo) have 
started [P3"| . If supplemented by a sensible continuation into the higher-twist 
region, a z distribution is obtained in agreement with the HERA data, provided 
the above-mentioned initial-state problem is dealt with as well. The required 
colour-octet MEs need not be much smaller than the naive estimates in table [I] 
IB3I. 



8 Conclusions 

The NRQCD factorization approach gives quarkonium production cross sections 
as a (finite) sum of short- distance coefficients times long-distance MEs. For the 
factorization to hold, the latter have to be process-independent. In this paper 
we have tested the universality of charmonium MEs by comparing their values as 
extracted from various reactions. Particular emphasis has been put to elucidate 
the uncertainties in such determinations. The results can be summarized as 
follows. 

Up to now, the cleanest evidence for the need of short-distance production 
of colour-octet states in the formation of 1 mesons comes from the J/ip and 
ip' p± distribution at the Tevatron. The numerical values of the corresponding 
colour-octet MEs are in line with the expectation from velocity scaling. Our 
investigation suggests that (Of ( 3 Si)) is somewhat larger and A4f_ 5 considerably 
smaller than currently extracted values (table [|; ip denotes J ftp or if)'). 

In the case of J ++ mesons: a colour-octet ME, (Cg c0 ( 3 S'i)), enters already at 
leading order in v 2 . Its numerical value as extracted from the Tevatron is highly 
correlated with the value of the colour-singlet ME (Oi cl ( 3 P±)) , whose value is 
less well known from potentail-model calculations than the colour-singlet if) ME. 
Indeed, B-meson decays into Xcj suggest a factor-of-2 smaller (Og c0 ( 3 Si)) value 
and a factor of 3 larger value for (C^ 1 ( 3 -Pi)). 

Although the full calculation of T decays into charmonia does not yet exist, 
the decay T — > J/ip + X does already provide upper limits on the J/ip colour- 

6 Here z is defined analogously to z, i.e. z (z) is the photon-energy fraction carried by the cc 
pair (the J/^>) in the proton rest frame. 
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octet ME. If the colour-singlet contributions were really as small as suggested 
by their a s scaling, this would then not only prove that short- distance colour- 
octet production is at work but also result in a (O^ ( 3 Si)} value larger than the 
current Tevatron estimate given in table |l[ 

Similarly, Z° data at LEP are about a factor of 3 above the prediction 
based on colour- singlet production alone. Inclusion of colour-octet processes 
with (0% ( 3 Si)} as currently deduced from the Tevatron reconciles nicely the 
agreement with the LEP data. 

Confronting theory and experiment of electron-positron annihilation into 
J / ipX at two different low energies is in principle the best way to separate the two 
colour-octet MEs (O^* (} S )) and (0 J 8 /lp ( 3 P )) that enter, in practically all other 
reactions, only in the linear combination (||) with almost constant a. However, 
current data suffer from statistics that are so low that they are compatible with 
pure colour- singlet production within 2 standard deviations. Neglecting the ex- 
perimental errors, a value of .M^ef is found somewhat smaller than the Tevatron 
one (table p. 

More accurate determinations of NRQCD MEs from e + e~ annihilation and Z° 
decays are mainly hindered by low statistics, and improvements can be expected 
with more data in the future. In contrast, constraints from T decays are currently 
limited by theory. The situation is similar for B decays into J /if} and if)': The 
largest uncertainty arises from the Wilson coefficient 2C + — C_ entering the 
effective weak Hamiltonian, which even in next-to-leading order is theoretically 
known only poorly. Pushing it to its maximum value, if> production is compatible 
with pure short- distance colour- singlet production. For most of the uncertainty 
range, however, the current Tevatron estimates for (0$ ( 3 Si)) are well compatible 
with data. On the other hand, 5 is once again found to be smaller than the 
Tevatron value in table [I], by about a factor of 2. 

Finally, photoproduction of J /if) as measured at HERA does not pose a prob- 
lem to NRQCD factorization, once colour-octet MEs compatible with all current 
data are being used, the p± broadening due to initial-state gluon radiation and 
intrinsic k± is taken into account, and the leading-order v 2 calculation is improved 
with the help of shape functions in order to extend its validity beyond a value of 
z of about 0.7. 
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